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Abstract 

Self-consistent hydro-dynamic one-loop quantum corrections to the Gross-Pitaev- 
skii equation due to the interaction of the condensate with collective excitations are 
calculated. It is done by making use a formalism of effective action and ^-function 
regularization for a contribution from Bogoliubov particles in hydrodynamic approx- 
imation. It turns out to be possible to reduce the problem to the investigation of a 
determinant of Laplace operator on curved space where a metric is defined by den- 
sity and velocity of the condensate. Standard methods of quantum gravity let us get 
the leading logarithmic contribution of the determinant and corresponding quantum 
corrections. They describe an additional quantum pressure in the condensate, local 
heating-cooling and evaporation-condensation effects of the Bose-condensed fraction. 
The effects of these corrections are studied for the correction from excited states in 
equilibrium situation. Response functions and form factors are discussed in the same 
approximation . 

1 Introduction 

The recent success of the experimental observation of Bose-Einstein condensation for sys- 
tems of spin polarized magnetically trapped atomic gases at ultra-low temperatures |], 0, ||] 
and futher investigation of their collective properties @] are stimulating development of the 
theory of a nonuniform Bose condensate and its collective excitations. Moreover, new excit- 
ing problems such as a description of the evolution of Bose condensate from relaxed trap, 
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dynamics of a collapse of the condensate for Li 7 atoms, heating-cooling phenomena, vari- 
ous coherence effects for the condensate and so on are being posed both theoretically and 
experimentally. 

A lot of theoretical efforts were spent to solve these problems. Evolution of the condensate 
from relaxed trap was considered in a number of papers both for symmetric and asymmetric 
cases P, ||, 0- All of these calculations were performed in Hartree approximation, i.e. based 
on Nonlinear Schrodinger equation or Ginzburg-Gross-Pitaevskii equation (GPE) [§]. It 
is obvious that the GPE describes strongly interacting Bose condensed gas well since it is 
possible to ignore non-condensed fraction. However next three points should be considered: 

1. At least for experiments with Rb 87 number of atoms is not such that strong interaction 
approximation (or Thomas-Fermi) approximation is sufficiently accurate to produce a 
robust ground state profile and realistic values of experimental observables. 

2. When the condensate expands after the release from a trap the density inevitably 
decreases and quantum corrections due to the interaction with particles in excited 
states should be considered even if the original density was high enough to be described 
by GPE. 

3. Effects like heating or cooling of the system and the collapse of atoms with attraction 
such as Li 7 (TTJ, [11, [12|] cannot be described even in principle by GPE. 



Therefore excited states and their influence on the condensate require consideration. 

Firstly, excitations of a trapped system, using the technique of Bogoliubov transforma- 
tion, have been considered by Fetter |13| (see also [14|). The spectrum and wave functions 



of excited states were derived by Stringari |T3| in the hydrodynamic approximation where 
the low energy excited states are sound waves in inhomogeneous condensate media. An 
analogous spectrum was obtained by Burnett and coauthors [TBI who investigated response 



of the condensate to a time-dependent perturbation. Using Stringari's approach, Wu and 
Griffin [[H] proceeded by quantizing the low energy excitations in a trap and diagonalizing 
the hydrodynamic Hamiltonian in terms of normal modes associated with the amplitude 
and phase. The depletion of the condensate and observable quantities (such as response 
functions) are some series over all excited states. It means that they cannot actually be 
evaluated analytically since contributions of various types of excited states are the same 
order of magnitude (as can be seen from Fig.l of Ref . ||1 7\\ ; see also discussion in section 5). 
Hence an effective method of summation as well as the generalization of the treatment to 
time-dependent and untrapped cases are needed. These problems are solved in the present 
paper (see also |18 



In recent paper Ref. [[[Jl the evolution of a Bose condensed gas under temporal variation 
of the confining potential was considered. It was found that there exist exact scaling relations 
for a 2-dimensional gas trapped in a harmonic potential whose frequency depends on time. 
Approximate relations also exist in the 3-dimensional case if the interaction is very small 
or very large (such that the kinetic energy can be neglected). This allows the authors to 
discuss principal questions of the condensate evolution from a relaxing trap and the depletion 
of the condensate. However this elegant work is restricted by the requirement of a harmonic 
potential confining potential. 
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We want to note also work [20| which considered nucleation of the condensate. It is 
possible to say that the present work is devoted to the consideration of opposite case where 
the condensate may evaporate. 

In the paper we derive a generalization of Gross-Pitaevskii equation for the condensate 
fraction taking into account one loop quantum corrections caused by the interaction of the 
condensate with the noncondensed fraction. However, there is other possibility to generalize 
GP equation. Indeed, GP equation in mean field approximation gives dynaimcal equation 
for both condensed fraction and full density of the gas since the deplition is ignored in the 
approximation. One loop correction for full density is obtained in Ref. JL8J that together with 
the results of the present paper allows to calculate the deplition in one loop approximation 
simply as a difference of the full density and density of the condensate fraction. 

The paper is organized as follows. In next section we consider the effective action and de- 
duce the most significant terms. In Section 3 we introduce functional integration to calculate 
the effective hydrodynamic action for the system and define one-loop quantum corrections 
which are responsible for the interaction of the condensate with excited modes. These quan- 
tum corrections are evaluated using ^-function regularization for functional determinants 
and an explicit and very simple expression for the one-loop quantum corrections is pre- 
sented. This permits in Section 4 the derivation of the quantum corrections to the equations 
of motion for the condensate in closed and explicit form readily used in numerical calcu- 
lations. Here corrections to the Nonlinear Schrodinger equation (due to the interaction of 
condensate with non-condensed fraction) are obtained. Using interpolation between low and 
high density limits, we show that the contribution from excitations are of order of the main 
terms in a region around the classical turning points and can define the condensate profile 
there. Section 5 is devoted to a consideration of Green functions and response functions. 
It is shown that all Green functions in the 1-loop approximation can be obtained from the 
effective action of Section 3 by functional differentiation. We conclude the paper with a 
discussion of the possible applications of our results. Two appendices contain basic facts 
and some technical details of ^-function regularizations of determinants of elliptic operators 
and methods of covariant variational differential calculus. 

Let us emphasize that our analysis does not depend on the details of the confining poten- 
tial and can be used in a variety of problems. The only simplification is usual hydrodynamic 
approximation. 



2 Exact dynamical equations for the condensate — re- 



view 



In this section we discuss relevant contributions from the non-condensed fraction which 
have to include in the exact equations of motion for the condensate. The estimate of the 
contributions is more transparent if we follow the usual approach of functional integral for 



the condensate dynamics [21 |. 

The Hamiltonian of the system can be written as 

2m 



H 



dV 



V^> + W> + v(x, t)4> + ij + iTii^^ipij 
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where v(x, t) is an external potential. For the case of trapped atoms the potential is supposed 
to be harmonic and can be varied or even switched off at some moment, say t = 0. More 
precisely it means that for the case of magnetically trapped atoms v = 9(—t)aT L (t)(p 2 + 
\ 2 (t)z 2 ), p 2 = x 2 +y 2 , (a± is oscillator lengths) has to be put in Eqn([J). Another parameter 
m Eqn(H) is I — the s-wave scattering length for atoms in the system p^| . The operators 
are Bose creation and annihilation operators. The corresponding vacuum- vacuum 
transition amplitude is given by the following functional integral: 



Z(p)= /D^ + (x,r)D^(x,r) exp[iS/h] 



with the action 



S 



dr dV { ihij)' 



2m 



2,771 

Vip + Vip + {v Tp)ip + ip + 4:irlijj + if} + ijjil> 

n 



(2) 



(3) 



where the parameter p is the chemical potential, controlling the number of particles in the 
system. 

As it is well-known the appearance of the Bose condensate is equivalent to the fact that 
besides "normal" quantities some "anomalous" ones have nonzero values. The principle 
"anomalous" quantities are one-point Green's functions (ijj(t)) and (ip + (t)) which are the 
amplitude of the condensate and have to be taken into account accurately. The standard 
approach to deal with the condensate is to make a shift of the fields ip + ,ip by some func- 
tions a + , a choosing from the condition of the vanishing of the "anomalous" Green's func- 
tions (vanishing of the one-point function provides analogies vanishing of other "anomalous" 
Green's functions). Formally, we consider the change of variables 

ip + = a + + tp + , ip = a + (p , 

which leads to a new form of the action (|3|): 

S = S 1 + S 2 + S 3 + S 4 + S 5 + S 6 (4) 



Si 
S 2 
S 3 
S 4 

s 5 



ida h 



2 r 



dr dV { ha 
00 J or 2m 



—00 
00 

-00 
00 



dr J&vlhf 



dr 2m 
ida h 2 



[ dr 2m 

dr [ d v \ha^-f 
00 J or 2m 



r { h 2 

dr dV{ -4ttZ 

00 J { 2m 

30 r ( h 2 
dr / dV { -Airl — 
00 J 2m 



Va + ■ Va + (v - 



^ 2 p)a a + ATcl(a a 



■ V<p + (v - —pii)^^ 



Vy? + • Va + (u - 



2m 
2m 
2m 



(5) 



h 

2m 



2 p)(f < a + 8nl(f < a aa 



Va + • Vip + (v 2~/ i ) a+ V 9 + 8nla + a + aip 



Aa + aip + (p + aatp + (p + + a + a + (ptp 



(Lp + (p) 2 + 2aLp + (p + <p + 2a + Lp + (p(p 



(6) 
(7) 
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By the definition of the shifts a + , a "normal" quantities for (p + ,ip — fields have to be zero, 
i.e. 



(ip) = J Dip + (x,r)D(p(x,T) {pexp[iS/h] = 

and 

(ip + ) = /D^ + (x,r)Dv?(x,r) (p + exp{iS/h] = . 



The last two equations are the main dynamical equations for the condensate wave functions 
(shift functions): 

7 + T = , 7 + f = (8) 

where 

., da h 2 _ 2 . 2m , , . . 

7 = —in— V a + (v y-u)a + Sirla aa (9) 

or 2m h 

and T is the sum of all diagrams with only one external outgoing line. Equation (|8]), in the 
mean-field approximation (neglecting T), becomes the Gross-Pitaevskii equation since the 
form of 7 (Q) is: 

da Ti 2 2tti 

— ih— V 2 a + (v n-u)a + 8nla + aa = . (10) 

or 2m h 

In previous work on Bose condensates in traps only equation (p~Q|) has been studied. 

Let us now estimate corrections to T from the interaction of the condensate with the non- 
condensed fraction. For the case of magnetically trapped atoms, i.e. a harmonic potential 
v = aj 4 (p 2 + \ 2 z 2 ) ( A is the anisotropy parameter) there is natural unit of length in the 
system, which is a±. After rescaling the coordinates in units a±, the interaction constant Anl 
is changed to g = 4irl/a± (in experimental conditions g ~ 0.05). Starting with the mean field 
solution for the condensate fraction in the trap, we can estimate effective interactions (Q) 
and (0). Indeed, the mean field solution for a number of particles iV ~ 5000 at the center of 
the trap is of order of 10. From this it follows that at the initial stage, the main contribution 
to the Gross-Pitaevskii equation is the term 2g\a\ 2 a ~ 100. However the next contribution 
from the interaction with non-condensed excited states with vertex 4g|a| 2 ~ 20 ^> 1 cannot 
be neglected. By comparison all other terms are quite small (g\a\ ~ 0.5, g ~ 0.05) and 
correspond to the weakly interacting non-condensed Bose gas contributions (for the original 
experimental setting with N ~ 2000 the numbers are 2g|a| 2 a ~ 30, 4g|a| 2 ~ 10, g\a\ ~ 0.3, 
g ~ 0.05). Hence the description of the transition stage from the original trapped condensate 
to the dilute phase require taking into account interaction effects due to terms (|) while 
the interactions (^) can be ignored. Moreover, cooling-heating effects, dissipation from the 
condensate and so on come from the interaction terms @ and cannot be neglected in the 
resulting description of the dynamics of the condensate. 

It is easy to see that the interaction (H) together with the quadratic part @ correspond 
to a Gaussian approximation in the integration over oscillations around mean field solutions, 
i.e. Bogoliubov's excitations. It is quite hard to integrate over Bogoliubov's excitations 
in general but this significantly simplifies in the hydrodynamic limit. As it was shown in 



T7H it is possible to obtain equations for hydrodynamic excitations straightforward from 



Bogoliubov's equations and then second quantize the excitations. Instead we will derive the 
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hydrodynamic description directly from the functional integral formalism, which simplifies 
the description. 



3 Quantum corrections to mean field equation in hy- 
drodynamic regime 

We are looking for the effective action which describes all physical quantities for the system 
(for easy introduction to the formalism of the effective action see, for example, Appendix 
A of Ref.[l8|]). For example, the effective action provides all Green's functions in the same 
approximation used to calculate the effective action itself. We will obtain it in a hydrody- 
namic one-loop approximation. To clarify the description hydrodynamical variables, density 
and velocity, should be used. More precisely, we start with the action (||) 



S 

h 



dr JdVi ihip' 



dip _ h 2 
dr 2m 



2vn 

Vip + Vip + (v t^ + ^P + 4ttZ(^ + ^) 2 

n 



and rescale variables as r 



— r ti 



2m' 



p, I — > j- to get the following action form: 



S 

n 



dr dV { iip~ 



dip 
dr 



Vip + Vip + {v- p)ip + ip + l(ip + ipf 



'12) 



(We will measure the action in terms of h everywhere below). Now we change field variables 
to the hydrodynamic ones: 



iP(x, r) = <Jp(x, r)e***rt , iP + (x, r) = yj p(x, r)e^ 
Then the action (KT ) takes the form (up to a complete derivative term): 



S 



dr I dV | ?gp - ( V Jpf - p{V V f -(v- p)p - lp 2 



(13) 



(14) 



It is not difficult to see that the classical equations of motion for the action lead to the 
equations: 



(Vif) 2 - v + n - lip + — V 2 ^p = 
dr Jp 

_^ + 2V(V^-p) = 0, 



(15) 



(16) 

which are equivalent to the Gross-Pitaevskii equation |15| after the introduction of the 
velocity variable c = — 2V<£> instead of (p. In velocity-density variables these equations 
look as hydrodynamic equations for an irrotational compressible fluid: 



dc 



- + 2v-2p + Alp-—V 2 ^p^ 







(17) 
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g + V(cp) = 0. 



Now let us shift our variables in (|T2| ) by the zero-order mean field solution to find one 
loop quantum correction for the effective action: 



i> i> + x 



— > %p + + x 



and we will keep up only terms including the square of x an d X + terms. Then the action 
(|l^) transforms to S(ip, + Si(ip, X-, X + ) where the action S\ has the form: 



clr / dV |z X + ^ - [Vx + V X + {0 - (V<p) 2 + 2lp} X + X 



-l^PX + X + ~ le~ 2iv 



PXX 



(19) 



where hydrodynamic representation ([13]) were used and potential v was expressed through 
mean field equation of motion neglecting kinetic term. Hamiltonian corresponding to Si is 



H= dV 



Vx + Vx + W - (V^) 2 + 21 P }X + X - l^PX + X + - le'^PXX 



-2i<p 



(20) 



Now we fulfill the following canonical transformation: 



X = y/pe 



'V 



a 
2~p 



+ ia 



x + = Vp? 



up 



a 
2~p 



ia 



The hamiltonian takes the form 



H 



dV 



p(Vd) 2 + p(pa 2 + {l + ^}a 2 + o-(V^Vd) - a(V(fV&) 



1 

2 £ 



+ : -Aipaa - ^Aipaa + v ' r ' a 2 + f (Vff) 2 + 



4p- 

(Vp)^ 2 + 1 
4p' 



(VP) 2 

4p z ~ r " ' L6p 3 



raVpVa + 



4p 



a 2 + dVpVd 



VpV tp 



(21) 



Last eight terms vanish in the hydrodynamical limit since p is a large variable (V<p ~ ^JTp 
and (p ~ Ip). So we have 



H= dV 



p(Vd) 2 + pfia 2 + {l + — }<r 2 + ff(V^Va) - d(V<pVa) 

4p 



f22l 



The corresponding action (in feynmann functional integral) is 
Si = 



p(Va) 2 + p<pa 2 + {/ + ^}ct 2 + 2cx(Vy?Va) 



4p- 



(23) 
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Integrating in functional integral over the field a we get the action for a field only: 



Si = J°° dr J dV j-p(Va) 2 + /(^ - 2Vy?Va) 2 - py?a 2 



(24) 



where 



./ 



4/ + <p/p 

The corresponding determinant does not contribute. 

Now we will use the large parameter p. Let us introduce the following large number 



Po = max {p} such that p 



Hence 



PO 



1 and new space-time variables t = 4Zp T, y% = V^Po^i- 



Si = 



=^= f°° dt [dV \A^d„ad v a + p<pa 

64/W J-oo J 1 



(25) 



where 



v = 2Vy? , f = Alf = 



1 + 



and the matrix A has a form 



fpi P-M -fvy 2 -fyivs 
fy 2 -fyiv 2 p-fvl -fv 2 y 3 

V f v 3 -M^3 -/u2«3 P-fvl j 

with its determinant det A = — p 3 /. 

Our next step is to cast the action in the covariant form. To do this we introduce auxiliary 
metric such that 
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fj,l/ 



V-det (11^11) ' 

One can easy to find the covariant metric from the equation above: 



-det ( || Av ||) 



that leads to the form: 



1^1 



1/2 ^-3/2 



r L/2 p 



( -f fvi fy 2 fy 3 \ 

fyi P-J^l -fvy 2 -fpiva 

fy 2 -fyiv 2 p-fvl -fv 2 y 3 

V /«3 -fviv 3 -fv 2 v 3 p-fvj j 
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Finally we get the metric 



P~ 1/2 / 1/2 





-PIS + v 2 




V2 


^3 \ 






1 










v 2 





1 





V 


V3 








1 / 



with the determinant g = det (Hs^H) = —p ? 'f- 

In this metric the action takes a covariant form 



Si 



1 



1 



dy \ -g 



where 



dy 



E = 



got 



pep 



aOa + Ea 
1 



-9 



gd v + E 



a 



(26) 
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This is the central equation of the paper. 

Now the effective action is I\ = — |tr ln[(64p ^ 3 )' 1 ^ 2 ( i::i + E)} with Laplace operator 
□ = —-^jd^g^^—gdj,. Let us again make use the existence of large parameter in the 

system. Indeed, the multiplier (64p ^ 3 )~ 1//2 gives a possibility to express main contributions 
to the determinant such as (see Appendix A): 

tr ln[(64p / 3 r 1/2 (D + E)\ = tr ln[D + E] - ^tr ln(64p Z 3 ) ($o(° + E) - L(D + E)) 

~ --tr ln(64p Z 3 ) ($(,(□ + E) - L(U + E)\ 

since in our regularization tr lnp + E] is order of the zeroth Seeley coefficient &o{^ + E) = 
jdydr ^—g^Q. We will consider only cases where the number of zero-modes does not 
change. So in what follows we will drop the term with £(□ + E). Returning to the initial 
variables all curvature tensors and the metric are written in the "physical" variables (i.e. 
without tildes). Moreover ln(64p / 3 ) 3> ln(p) so that we can substitute ln(64p/ 3 ) instead of 
ln(64p / 3 ). Summarizing, we obtain an expression for the first quantum correction to the 
effective action: 

Ti = ^J dxdr v^o;in(64pZ 3 )^ (n + E) , (27) 



with the metric 



,V2 



/V 2 



and 



Vi = 2di<f , / 



-p/f + v 2 

Vl 
V2 
V3 

1 

4/ + (pjp ' 



Vl 

1 





V2 



1 





^3 \ 




1 



E 
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4 Quantum corrections to equations of motion 



In this section we derive from required quantum corrections to equations of motion 
(see, for example, |18|). To make the consideration self-consistent we should not differentiate 
ln(64p/ 3 ) because in our approximation it behaves like a constant. In that case, the action 
we have to vary is covariant again. We find: 



dg, u 6% | dE5% 
dp Sg^ u dp 5E 



5 J± 



iln(64p/ 3 ) Jd 4 y 



Sg^x) 5$ 5E(x) 5% 



5<p(y) Sg^(x) 5cp(y) 5E(x) 



Let us remember that 



(An 



— V 2 R + —R 2 - —R^R^ 
30 72 180 M 



— R, luap R^ ap + -RE + -E 2 - -V 2 E 
180 ' p 6 2 6 



We note that two of the terms above are irrelevant: 



| dx v^V 2 F = |dx ^j-^d^^jg^dvF = Jdx d^^jd^F) 



since this is just a complete derivative. Hence for our purposes it is sufficient to consider 
only 



$ (n + E) 



(An 



72 



Rm.R^ H R, n , rrn R p,l ' C!f> H — H — _£/ 2 

180 M 180 6 2 



Long but quite straightforward calculations lead us to the following result for the func- 
tional derivative (see Appendix B for details of the calculation): 



5g 



-$r 



(An) 



11 i 11 

-—R 2 R rTn R Tf} H R rTCin RR ul>OL ^ H — RE H — E 2 

72 180 p 180 paP 6 2 



+ q R ^uE + —R^R - — R^R'J* + —R papa R v ' Tpa - —g^ u UR 



+ 



1 



180 

taking into account 



— V p V a R 
45 



— $ (n + £) = -^V P T^ $ o(n + E) . 



5^ 
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This means that the equations of motions with 1-loop quantum correction take the form: 



v + (i - 21 p + -ln(64p/ 3 ) 



+ 



dp 5guv dp 5E 



^ + 2V(V^-p) + ^ln(64pZ 3 ) 



5% 



'dE5% 
dip 5E 



(2* 



, (29) 



These two equations (|28f) and (P9|) replace the mean field hydrodynamic equations ([lqJTq ) 
for the density and velocity of the condensate and are the main result of the paper. They 
contain contributions of the excited states to the dynamics of the condensed fraction and 
provide a description of effects such as depletion and heating-cooling. For example, we 
note that the additional quantum pressure is due to dependence of the Bogoliubov particles 
determinant on the density of the condensate while an evaporation of the condensate (see 
continuity equation (|29|)) comes from phase- dependence of the determinant. In particular, for 
the stationary situation when the equilibrium condensate velocity is equal to zero there is no 
evaporation of the condensate and the continuity equation holds for the fraction. This case 
was considered in Ref.fllBl 



However we think that equations ( Pq) and (|29|) are of the main 
interest in non-equilibrium problems such as evolution of the condensate under changing of 
trap shapes, collapse of condensate cloud and corresponding heating and evaporation of the 
condensate and so on. We hope to consider these problems in forthcoming papers. 



5 Green functions, response functions and the like 

Equations of motion and quantum corrections for them are the central questions of previous 
consideration. However in many applications related to experiment it is very important to 
analyze various response functions and form-factors. They can be expressed in terms of 
some combinations of Green functions of the theory. That is why in this section we shortly 
consider a calculation of Green functions in the effective action formalism. 



As it is shown in Ref [rj| the effective action (in any self-consistent approximation) allows 
the evaluation of all Green functions in the same approximation by just taking of variational 
derivatives. Since the effective action in 1-loop approximation was obtained the problem of 
Green function calculation is the problem of functional differentiation only |24 . 

We now formalize all said above and give formulas to evaluate Green functions in effective 
action formalism. Let W(J) being generating functional of connected Green functions. Then 
quantities 

ly „( Xl ,..., lJ = _i_..._i_ lyW (30) 

are connected Green functions in the external field J. The effective action is defined then 
by the Legendre transformation of W: 

r(a) = W(J(a)) - a J (a) , a = (31) 

oJ 
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where the function a(x) = (<p(x)) = W±(x] J) is the first connected Green function (p0|). It 
is easy to see that the functions a and A are related by the second relation (0). Then the 
quantities 

5 5 

r » (xl '-"' X " ) = ^"-*aW r(a) (32) 

define 1PI (one-particle-irreducible) Green functions for the theory with mean field a(x) = 
(tp(x)). Knowledge of 1PI Green function is equivalent to knowledge of any (whole or con- 
nected) Green functions for the corresponding system. For the first 1PI Green function we 
have from fl3~l|) 

r x (x) = -J{x) . (33) 

All connected Green functions can be expressed in terms of 1PI ones. Indeed, differentiating 
(j33|) J we obtain 

W 2 T 2 = -1 , W 2 = -r 2 1 (34) 

or, in expanded form, 

r dz W 2 (x, z)T 2 (z, y) = -5(x - y) . (35) 



Differentiating now the second relation in (|34[) on J and using the following rule 

5J 2 5a 

we can derive expressions for all the higher connected Green functions through the 1PI ones. 
For example, for the third connected function we have 



1 3 

1 2 



r 



or 



W 3 (xi,x 2 ,x 3 ) = - J dy!dy 2 dy 3 T 2 1 (xi,y 1 )T 2 l (x 2 ,y 2 )T 2 1 (x 3 ,y 3 )T 3 (y 1 ,y 2 ,y 3 ) (36) 

and so on. It means that since we know 1PI functions we have to solve the differential 
equation (j35D for the two-point correlation function and then find all other connected Green 
functions by integration of 1PI functions with two-point correlators as in Eqn. (|36|) . 

Putting J = in ([33]) we get equations of motion for ip and ip + (|28|j2^) . To obtain, for 
example, the Green function ((ip + (x) — (tp + (x)))(ifj(y) — (ip{y)))) we twice differentiate the 
effective action on substitute solution of equations of motion and finally invert the 

result in the sense of the kernel of a integral operator. 

As an example of usage of this technique we calculate two-point Green function in the 
mean field approximation (tree or 0-loop approximation). In this approximation the effective 
action has the form: 

T = y°° dr J dV l-i^r^ - W> + V?/> - {v - fM)i) + ip - /(^ + ^) 2 | . (37) 
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such that, in the equilibrium hydrodynamic picture, the matrix of second variational deriva- 
tives of the effective action can be written as 



/ 



5 2 T 



5 2 T 



6 2 T 



\ 



( 

V 



-t-^-V 2 -(v-p)-2l(^ + ^) 



5 2 T 6 2 T 

\ 5^(x)S^+(y) 8ip(x)5ip(y) ) 

. d 



" Tr ~ V 



2 (v-fi)- 2l(ip + ip) 
2l^ + ) 2 



where ip^ = ('tp^). The matrix of two-point connected correlation functions of fluctuations 

ijj(x) = tp(x) — (ip(x)) and ^ + (x) = — (^ + (x)) 



G 



(^(x)^(y)) ^ + {x)i,{y)) 
(tp(x)^+(y)} (tp(x)^(y)} 
is defined solution of the equation 



5 2 T-G 



-I . 



Solving the equation we get the following expressions for the correlators. One can check that 
the Green functions calculated in this manner coincide with the Green functions (p(x)p(y)), 
(p(x)ip(y)), (ip(x)(p(y)) analyzed in spectral representation by Wu and Griffin |I7[]. To include 
the first quantum correction (1-loop quantum correction) term (|27|) should be added to action 



In summary, various response functions, form-factors and Green functions are evaluated 
from effective action in any given approximation taking variational derivatives and integrat- 
ing them with the two-point correlation function. This can be done in general using formulas 
of this section. 



6 Conclusion 

In the paper a self-consistent approach to the calculation of 1-loop quantum corrections 
to the Gross-Pitaevskii equation due to the interaction of the condensate with collective 
excitations is considered. To do this the hydrodynamic approximation was used. In this 
approximation excitations are equivalent to sound waves on a condensate background. This 
opens the possibility to use methods of quantum gravity and theory of quantum gauge fields 
where the problem of calculation of effective action for gravitational and gauge backgrounds 
(or more precise, quantum corrections due to other quantum fields) is common problem. 

Many of the methods to approach the problem were developed during the last four 
decades. They are (^-function regularization for determinants of operators, Schwinger-De- 
Witt-Seeley expansion of heat kernels P7L E3, 123], covariant methods of calculation of Seeley's 



13 



coefficients and the covariant perturbation technique for effective actions pUL ffT]|. Although 



the methods are well-known in field theory they are not so familiar in condensed matter 
physics and the theory of coherent systems. An accurate account of corrections is very 
complicated even for a few first orders and the solutions are obtained by using the covariant 
perturbation technique and curvature expansions. On the way effects nonlocal in space 
and time appeared. In quantum gravity framework such effects play significant role in the 
gravitational collapse problem, Hawking radiation and so on. 

However to benefit from it sometimes very general and only basic information about 
principle scales in systems in question is required. Indeed, this is the only information needed 
to extract leading logarithmic contributions while the calculation of other corrections is much 
more complicated. That is why we think that the developed approach can be widely used 
in many problems which have not much to do with Bose-condensation of trapped atoms or 
liquid Helium. 

At the end let us stop on other applications of the presented method. It is easy to imagine 
other condensed matter examples of problems treatable by the same technique. Clearest and 
simplest of them are antiferromagnets, Josephson arrays and superconductors in nonuniform 
(in space and time) external magnetic fields. However many of other physical systems in a 
non-equilibrium background are potential field of applications. 
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Appendices 



A Determinants of elliptic operators 

Usually the first quantum correction to the effective action is ill defined. The point is it is 
divergent. This is just the well-known ultraviolet divergence of the quantum field theory. 
Indeed, we can rewrite the functional determinant as 

lndet A = Tr InA = ln[]A n = J]lnA n (A.l) 

n n 

where A n are the eigenvalues of the operator A. This series is easy to show to be divergent. 

That means that one needs a regularization. This point was investigated very thoroughly 
by many authors and it is found that in quantum gravity and gauge theories the most 
appropriate regularizations are the analytical ones. The functional determinants can be well 
defined in terms of the so called ^-function. 

At first let us use a Wick rotation to produce a Laplace operator from the wave operator. 
It allows to make use of methods for the evaluation of the determinants of elliptic operators 



| 25fl . Now we are ready to introduce the ^-function of an elliptic operator A: 

t(s,A) = J2K s 

i 

where {Aj} are eigenvalues of the operator A. Then 

Tr lnA = - ^-((s,A) 
To study the behavior of the ^-function it is common to use the formula 

oo 

KS = T(s)I dt tS ~ lex P(- tA ») • 

^ ' 

Summing over all nonzero eigenvalues of A we get 

oo 

C(s, A) = -L jdt t s -\Tr exp(-tA) - L(A)) (A.2) 



o 



where L(A) is a number of zero-modes of A. Representation ( A.2 ) is valid when the integral 



converges. For nonnegative operator it always does converge as t — > +00. Conditions of the 
convergency of the integral as t — > +0 depend on details of the operator A. For Laplace 
operator in 4D it converges as t — > +0 if Res > 2. 

There exists the well-known Seeley expansion for the Tr (exp(-tA)): 

Tr exp(-M) = £ §- k {A)t~ k + p(t) , 

fc>0 
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where \p(t)\ is bounded by a constant times t as t — > +0. are called Seeley coeffi- 

cients and play important role in the investigation of elliptic operators and their topological 
properties (for example in the Index theory). For example, for 4D Laplace operator in an 
external potential E = k > 3, and 

$_ 2 = (4tt)- 2 , 

$_i = — (4n)~ 2 (E + , where R — scalar curvature , 

and 

*° = ^ (-To V2R + T2 R2 - Wo R ^ + Wo R ^°" + l ER + ^ " \ V%E 

Here R^ vap and i?^ are Riemann and Ricci tensors correspondingly. In the paper we use 
the following definitions for them |2B] : 

and 

pA ^"^p/v ^"^PP I "per pA per pA 

P^P ~~ Q x p Q x u H v PO" W va i 

with the Christoffel symbols: 

rCT \ ua ( d 9nv dgx^ _ dg^x 
A " 2 9 1 &r A <9x" ax 1 - ' ' 
They are defined by a metric of a curved space. The latter coefficient $o is particular 
important for the calculation of determinants. 

Splitting (|A.2[ ) into an integral over [0, 1] and one over [1, +oo] we get 

« s > = rwl£— + — ; — + 



/CO rl 
dt Tr exp(-tA)^ 1 + J dt p(t)t s_1 



The singularity at s = turns out to be removable since lim s ^o s ^{s) = 1 and £(s) is defined 
by analytical continuation from the half-plane Re s > 0. 
From the relation lndet A = — C'(0) we find that 

$-fc(A) , , r/ ,,\ f°°dt r 



IndetA = ^_^U + r'(l)($ (A)-L(A))- f° Trexp(-M)- 
fe>0 K ' Ji t 

- f 1 ^ [Tr exp(-tA) - Y,^-k(A)t- k ) . 

This is (^-regularized In det A which we use in the paper. 

To conclude this section we give the relation between In det A and In det a A where a is 
a number parameter. It is easy to see that 

C(s,aA) = C(s,A) -a' 3 

because Xi(aA) = aXi(A). This leads to the relation: 

In det a A = In det A + In a ■ C(0, A) = In det A + In a ■ (® Q (A) - L(Afj (A.3) 

which we use intensively in the paper. 
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B Variational derivative of the effective action 



In this appendix we give details of the calculation omitted in section 3 to do not interrupt 
the main stream of the consideration. To make the calculation more efficient we make 
use covariant form of the quantum correction since there exists well-know way to simplify 



covariant variation calculation 26 



We are interested in the following variational derivative 

S 



$0 □ 



where 



(47T) 



dx a/^O — R 2 RnuR^ + Rma,, 

■ \72 180 M 180 M p 



R ^ap + _ RE + \ E 2 

6 2 



For variation of the volume element we have 

1 



where is defined by as 



-9 



Sg^g 



fw 



-gtiJg 



pi/ 



A' 



g-g 



and the following relation was used 



Sg^g uX 



For variations R 2 , R flu R' lu , R fluap R flucrp one can obtain 

5(RE) = 5R-E = 5g 12 R 12 E + g 12 5R l2 ■ E 
5R 2 = 25R-R = 25g 12 R l2 R + 2g 12 5R 12 ■ R 
6R 12 R 12 = S(g 13 g 24 R 12 R M ) = 25g 12 R 13 R 2 3 + 25R 12 R 12 
5Ri2saR 12M = 5(g 15 g 26 g 37 g 48 R 12U R 5678 ) =45g 12 R 1U5 R 2 345 + 26R 12U R 

where numbers 1, 2, 3, . . . mean indices fix, fi 2 , /X3, . . . respectively. Hence, 



1234 



8®r 



(47T) 



dx 



-g $g 



1 



1 



72 



R 2 D pl2 

180 



(B.l) 



(B.2) 



-R1234-R 1234 + —RE + —E z 



180 



6 



1 

45" 



1 
6 



(4tt 



—g 12 6R 12 R - —6R 12 R 12 + —6R 12U R 12M + -g 12 6R 12 E 
36 y 90 90 6 y 



So we need to calculate last three terms. 
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For variations of Ricci and Riemann tensors we have 
8R 



5R 



puap 



-g^5g 12 R 2vap + g^V p (5Tl a ) - g^V a {5Tl p ) 



(B.3) 



1 r. 



= -9/iiSg R 2vap + - [VpVJg^ + V p V v 5g pa - V p V p 5g v 
-V a V p 5g pu - V a V u 5g pp + V a V^Sg up 

Let us consider covariant divergence of some vector T k : 



(B.4) 



V P T» = V% = a M T" + T%T» 



1 



-9 



Hence 



is a pure divergence and can be dropped under integration. Using this we can easily derive 
the following expression for variations of the integrals we are interested in 



J dx ^gg^SR^R = J dx 



da; y^gSRuuR^ 



dx 



+-{V p ,V u }R-g pv UR 
-OR^ - g pu V 1 V 2 R 12 + ViV^i?, 1 + ViV.iT 1 



dx y/^SR^pR^" = dx V^gSg^ -R pl23 R u 123 + 2V 2 V 1 i? 



fil2v 



J dx y/^gg^SR^E = Jdx ^= 



Summing all terms we obtain 



--{W p ,V u }E-gp V UE 



6g 



-$r 



(4tt) 



-g\ ~2 9 ^ 



■—R 2 — R nn R ap H R rrnn RR' Tpa ^ H — RE -\ — E 2 

72 180 p 180 p ^ 6 2 



+ qR^E + —Rp U R - —R pa R u a + —Rp„p a R u apa - —g pu OR 
+ ^{Wp,W u }R-^g piJ nE + ^{Wp^ v }E 



+ 



1 

180 



-2V°V u R pu + UR pv + g pv V a V p R, 



up 



+ 45 ^V% P , 
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